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ABSTRACT 

Observations of low mean metallicity of damped Lyman-alpha (DLA) quasar absorbers at all red- 
shifts studied appear to contradict the predictions for the global mean interstellar metallicity in 
galaxies from cosmic chemical evolution models. On the other hand, a number of metal-rich sub-DLA 
systems have been identified recently, and the fraction of metal-rich sub-DLAs appears to be consid- 
erably larger than that of metal-rich DLAs, especially at z < 1.5. In view of this, here we investigate 
the evolution of metallicity in sub-DLAs. We find that the mean Zn metallicity of the observed sub- 
DLAs may be higher than that of the observed DLAs, especially at low redshifts, reaching a near-solar 
level at z < 1. This trend does not appear to be an artifact of sample selection, the use of Zn, the 
use of iVni-weighting, or observational sensitivity. While a bias against very low metallicity could 
be present in the sub-DLA sample in some situations, this cannot explain the difference between the 
DLA and sub-DLA metallicities at low z. The primary reason for the difference between the DLAs 
and sub-DLAs appears to be the dearth of metal-rich DLAs. We estimate the sub-DLA contribution 
to the total metal budget using measures of their metallicity and comoving gas density. These calcu- 
lations suggest that at z < 1, the contribution of sub-DLAs to the total metal budget may be several 
times that of DLAs. At higher redshifts also, there are indications that the sub-DLAs may contribute 
significantly to the cosmic metal budget. 
Subject headings: Quasars: absorption lines-ISM: abundances 



1. INTRODUCTION 

As the cycle of stellar birth and death progresses in 
galaxies, the interstellar gas is expected to become in- 
creasingly metal-rich. A fundamental prediction of most 
chemical evolution models is that the mass-weighted 
mean interstellar metallicity of galaxies should rise from 
a low value at z > 2 to a near-solar value at z = 
(e.g., Pei et al. 1999; Somerville et al. 2001; Cora et 
al. 2003). Indeed, the mass- weighted mean metallicity 
of local galaxies is close to solar (e.g., Kulkarni & Fall 
2002; Fukugita & Peebles 2004). How was this present 
level of metallicity reached? A direct way to study this 
enrichment history is through the absorption lines super- 
posed by galaxies on the spectra of background quasars, 
in particular the DLAs (defined to have neutral hydrogen 
column densities log Yhi > 20.3) and the sub-DLAs (de- 
fined to have 19.0 < log N m < 20.3; Peroux et al. 2003). 
These absorbers contain sufficient neutral gas mass to ac- 
count for a significant fraction of the visible stellar mass 
in present-day galaxies (e.g., Storrie-Lombardi & Wolfe 
2000; Peroux et al. 2003, 2005; Prochaska et al. 2005). 
Indeed, hydrodynamical simulations naturally produce 
DLAs from cool gas accumulating at centers of cold dark 
matter (CDM) halos and predict DLAs to be closely as- 
sociated with star forming regions (e.g., Haehnelt et al. 
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2000; Cen et al. 2003; Nagamine et al. 2004). 

Previous studies of element abundances in high-iVm 
quasar absorbers have focussed primarily on DLAs, es- 
pecially those at high redshift (e.g., Prochaska et al. 
2003b). Recently, we have started to expand the low- 
z DLA Zn data and tripled the z < 1.5 sample so far 
(Khare et al. 2004; Kulkarni et al. 2005; Meiring et al. 
2006a, 2006b; Peroux et al. 2006b). These studies find 
that the DLA mean metallicity stays substantially sub- 
solar even at low redshifts and increases relatively slowly 
(see also Prochaska et al. 2003a). Furthermore, based 
on deep emission-line imaging studies (e.g. Kulkarni et 
al. 2006 and references therein), the star formation rates 
(SFRs) of a large fraction of DLAs appear to fall far be- 
low the predictions from the cosmic star formation his- 
tory (inferred from galaxy imaging surveys such as the 
Hubble Deep Field, e.g., Madau et al. 1998). A missing 
metals problem has also been noted in high-z DLAs from 
SFR estimates based on C II* absorption (Wolfe et al. 
2003). These results could potentially present important 
challenges to our understanding of galaxy evolution. 

Why do DLA data disagree with the mean metallic- 
ity in galaxies predicted from the global star formation 
history? Smoothed particle hydrodynamics (SPH) sim- 
ulations (Nagamine et al. 2004) indicate that true DLA 
metallicities could be 1/3 solar at z = 2.5 and higher 
at lower z. Has the DLA metal content been underesti- 
mated, or is there a larger quantity of metals in a hitherto 
unexplored absorber population? While dust selection or 
metallicity gradients can lead to a flat metallicity-redshift 
relation, the extent of these effects is still uncertain (e.g., 
Boisse et al. 1998; Vladilo & Peroux 2005; Ellison et al. 
2005; Akerman et al. 2005; Chen et al. 2005; Zwaan 
et al. 2005). It has been suggested recently that dust 
obscuration is not significant and that metal-rich DLAs 
are inherently rare (Meiring et al. 2006a; Khare et al. 
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2006, hereafter K06; Herbert-Fort et al. 2006). In view 
of this, it is important to explore other possible sites for 
the missing metals. 

Over the past two years, we have been investigating 
the metallicities of sub-DLAs at z < 1.5 (Khare et al. 
2004; Peroux et al. 2006a, 2006b; Meiring et al. 2006b). 
These studies have led to the discovery of five solar or 
supersolar metallicity sub-DLAs at 0.7 < z < 1.5. Before 
these studies began, only a single metal-rich sub-DLA 
was identified (Pettini ct al 2000). Indeed, five of the 
nine sub-DLAs with Zn data at z < 1 and six of the 17 
sub-DLAs at z < 1.5 have solar or supersolar metallicity. 
Despite the small number of measurements (caused by 
the focus on DLAs in the past), the large fraction of 
metal-rich sub-DLAs suggests that log AThi ~ 20.0 may 
be the optimum range for revealing metal-rich systems, 
since such 7Vhi is high enough to have a significant metal 
column density and low enough to make dust attenuation 
less important (Lauroesch et al. 1996; Vladilo & Peroux 
2005; Peroux et al. 2005, 2006a). To investigate this 
further, here we examine the evolution of the observed 
amount of metals in neutral gas in sub-DLAs and their 
contribution to the cosmic metal budget. 

2. METALLICITY EVOLUTION OF DLAS AND SUB-DLAS 
2.1. Observed Trends 

We now compare the observed Zn metallicity evolution 
of DLAs and sub-DLAs. For several reasons, the nearly 
undeplcted element Zn is the most robust metallicity in- 
dicator in DLAs. Zn tracks Fe closely, to within ~ ±0.1 
dex, in Galactic halo and disk stars for [Fe/H] > —3 (e.g., 
Cayrelet al. 2004; Chen, Nissen, & Zhao 2004). Further- 
more, Zn II has two, often unsaturated, absorption lines 
at AA2026, 2062 that, for low-z QSOs, are free of blend- 
ing with the Lyman-a forest. 7 We do not include [O/H] 
estimates based on X-ray absorption measurements be- 
cause these may be systematically biased with respect to 
[Zn/H]. Finally, we do not include metallicity estimates 
based on emission lines detected from absorber galaxies 
(available only for a few DLAs) since these do not nec- 
essarily probe the metallicity of the absorbing gas. 

Our samples contain 119 DLAs with 0.1 < z < 3.9 
and 30 sub-DLAs with 0.6 < z < 3.2. (No Zn data exist 
for sub-DLAs with z < 0.55 because Zn II AA2026, 2062 
can be accessed at z < 0.55 only with space-based UV 
spectrographs, and the few existing UV measurements 
at z < 0.55 focussed on DLAs rather than sub-DLAs.) 
The samples are based on our recent data (Khare et al. 
2004; Kulkarni et al. 2005; Meiring et al. 2006a; 2006b; 
Peroux et al. 2006b) and data from the literature. Re- 
cently Prochaska et al. (2006) reported 2 metal-rich sub- 
DLAs at z ~ 1.8, but the use of these systems for study- 
ing the sub-DLA metallicity evolution could be debated 
since they were discovered in a study targeted specifi- 
cally toward metal-rich systems, i.e. those with strong 

7 Some authors have advocated using S and Si. We prefer to 
use Zn alone to have homogeneous data without the complication 
of nuclcosynthctic differences between different elements. Further- 
more, Si can be considerably depleted in some cases (e.g., [Si/Zn] 
= -0.9 dex in the cold diffuse interstellar clouds). Finally, the only 
accessible lines of Si and S are often saturated and the S lines nearly 
always lie in the Lyman-a forest. Nevertheless, as discussed below, 
we have checked that using S and Si measurements in cases of Zn 
limits makes essentially no change to our results. 



Zn lines detected in the SDSS spectra. We therefore do 
not include the Prochaska et al. (2006) systems in our 
primary sample, but do discuss the effect of adding these 
two systems. For the same reason, we do not include 
the DLAs from Herbert-Fort et al. (2006), but discuss 
the effect of adding these DLAs later. Table 1 lists our 
primary sub-DLA sample (sample 1), together with the 
DLA sample. The DLA sample contains 68 detections 
and 51 limits, while the sub-DLA sample contains 13 de- 
tections and 17 limits. 80% of the DLA Zn limits and 
53% of the sub-DLA Zn limits are at z > 1.5. 

Our analysis uses statistical techniques described in 
Kulkarni & Fall (2002). We divided the 119 DLAs in 
6 redshift bins with 19 or 20 systems each and the 30 
sub-DLAs in 3 redshift bins with 10 systems each. In 
each bin, we calculated the Nni-weighted mean metallic- 
ity, using the Kaplan-Mcicr estimator to account for the 
presence of Zn II upper limits. These calculations were 
done using the astronomical survival analysis software 
ASURV Revision 1.2 (Feigclson & Nelson 1985; Isobe & 
Feigelson 1990). Fig. 1 shows the results. The circles de- 
note DLAs, while the squares denote sub-DLAs (sample 
1). The look-back times are estimated assuming the con- 
cordance cosmology Q m = 0.3, Qa = 0.7, and H = 70 
km s _1 Mpc -1 . In each bin, the data point is plotted at 
the median look-back time, with horizontal bars indicat- 
ing the range of look-back times spanned by the systems 
in that bin. The vertical error bars denote the 1 a uncer- 
tainties in the Nju-weighted mean metallicity, including 
sampling as well as measurement uncertainties. 

Clearly, the global mean metallicity of DLAs is sub- 
stantially sub-solar even at low z (—0.82 ± 0.15 in the 
bin 0.1 < z < 1.2). Since the lowest redshift bin spans a 
large time interval, we also calculated the metallicity in 
two halves of this bin with 10 DLAs each and obtained 
nearly identical values (denoted by dashed circles in Fig. 
1). On the other hand, the two lower- z sub-DLA points 
appear higher than the DLA points. The Nui-weighted 
mean metallicity in the sub-DLA bin at 0.6 < z < 1.4 is 
0.05 ± 0.27. A linear regression fit to the Nni-weightcd 
mean metallicity vs. redshift data gives a slope of -0.25 
± 0.07 for DLAs (in agreement with Prochaska et al. 
2003a and Kulkarni et al. 2005), and -0.36 ± 0.19 for 
sub-DLAs. 

We next attempt to judge the dependence of our results 
on the sample definition, the use of Zn, and the use of 
iVm-weighted metallicities. We repeated our calculations 
for a number of changes in technique and assumptions. 
These results are summarized in Table 2. Concerning 
sample definitions, for example, we find that including 
the 2 metal-rich sub-DLAs reported by Prochaska et al. 
(2006) increases the metallicity in the middle redshift bin 
by 0.17 dex. We also find that the DLA results change 
very little if the DLAs from Herbert-Fort et al. (2006) 
are included: metallicities within 0.01-0.10 dex of those 
shown in Fig. 1 for 5 of the bins, within 0.2 dex for 
the remaining bin at 2.6 < z < 3.1, and the same lin- 
ear regression slope (-0.25 ± 0.07). If [S/H] or [Si/H] 
are used, where available (z > 1.8 and z > 1.3, respec- 
tively), in cases of Zn limits, and solar S/Zn or Si/Zn 
ratios assumed, the 2 higher redshift sub-DLA points in 
Fig. 1 change by < 0.05 dex. (There was no change in 
the lowest redshift bin where S or Si data do not exist, 
and the relative number of Zn limits is anyway lower.) 
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The resultant slope of the sub-DLA metallicity-redshift 
relation is virtually identical to that obtained from Zn 
alone, i.e., —0.39 ± 0.19. The sub-DLA trend suggested 
by Fig. 1 is, therefore, not an artifact caused by the 
use of Zn. Finally, repeating our calculations using the 
unweighted mean metallicity also gives similar results. 

It thus appears that the sub-DLA global mean metal- 
licity may be higher than that of DLAs, reaching a near- 
solar value at low z, consistent with predictions for the 
global mean metallicity in galaxies from chemical evo- 
lution models (see the curves in Fig. 1). Of course, it 
should be remembered that not every sub-DLA is metal- 
rich, and that our result holds for the Ajji-weightcd mean 
metallicity This result agrees with the conclusion of 
York et al. (2006) that the average [Zn/H] > -0.2 for 
a composite spectrum made from 698 SDSS absorbers 
with average log Nm ~ 20.0. Furthermore, our result 
supports the conclusion of K06, based on observed cor- 
relations between [Zn/H] and Nm, and a detailed con- 
sideration of various selection effects involved therein, 
that sub-DLAs are more metal-rich than DLAs. Our re- 
sult also supports a similar suggestion by Peroux et al. 
(2003) based on [Fe/H]. 

It is important to note that the low-z sub-DLAs in our 
sample were not chosen a priori to have strong Zn lines. 
[In fact, none of the 10 SDSS sub-DLAs with z < 1.5 in 
our sample showed Zn detections in the low-resolution 
SDSS spectra. Three of these 10 SDSS systems show 
detections of other relatively weak lines (Si II A 1808 or 
Fe II A2249) in the SDSS data, but interestingly, all of 
these systems have low metallicity [Zn/H] < — O.8.] The 
low-z sub-DLAs that we observed for Zn II were chosen 
based on a known H I column density from HST UV data, 
moderately strong Fe II A2344 lines in the SDSS spectra, 
and a reasonably bright background quasar (usually with 
V < 19 magnitude) with coordinates appropriate for the 
observing site and time. It is important to note that such 
systems are much more common than the "metal-strong" 
systems (with strong Zn II A2026 or Si II A1808 lines seen 
in SDSS spectra) used in the sample of Herbert-Fort et al. 
(2006). Furthermore, the samples in Fig. 1 include data 
from the literature as well (50% of sub-DLAs and 43% 
of DLAs at 0.6 < z < 1.5). Finally, the DLAs at z < 0.5 
from Kulkarni et al. (2005) were selected regardless of 
the Fe II strengths, and most were still found to be metal- 
poor. 

2.2. Some Potential Selection Effects 

2.2.1. Does Mg II selection bias the low-z sub-DLA sample 
toward metal-rich systems? 

Most of the low- z sub-DLAs and DLAs for which we 
have obtained Zn measurements were discovered in the 
HST UV surveys of Mg Il-selected systems (Wtfg > 0.3 
A; e.g., Rao, Turnshek, & Nestor 2006). One might, 
therefore, wonder whether the Mg Il-selected systems 
would tend to be more metal-rich simply because they 
are selected by the presence of strong Mg II lines. If a 
constant rest equivalent width threshold of Mg II A2796 
implied a fixed threshold of AM g n, then one may think 
that systems with low JVhi would be included only if they 
were more metal-rich, compared to systems with higher 
iVni- However, there are several indications that this 
does not appear to be the reason for the higher metal- 



licity we see in sub-DLAs compared to DLAs. First, 
Figs. 5-7 of York et al. (2006), based on hundreds of 
SDSS Mg II systems, show that strong W£fg$ occur for 
all E(B — V). This, together with the trend of decreas- 
ing metallicity with increasing E(B — V) in the data 
of York et al. (2006) (e.g., their Table 7) implies that 
using systems with large does not bias one to- 

ward selecting only metal-rich sub-DLAs. Second, the 
Mg II A2796 line is almost always highly saturated in 
DLAs and sub-DLAs. The systems with larger 
have larger velocity spreads rather than larger Mg II col- 
umn densities. For example, the compilation of Ledoux 
et al. (2006) for DLAs/sub-DLAs based on VLT UVES 
spectra shows that the velocity spread of such systems is 
typically about 100 km s _1 , and is even > 200 km s _1 
in some cases. In other words, selecting systems with 
larger W^gl does not necessarily translate into selecting 
systems with larger A^Mgii- Finally, as seen in Table 1, 
the JVhi values of the sub-DLAs in our samples are rela- 
tively large: 86% of the 30 systems in our primary sam- 
ple have log N m > 19.7 and 73% have log 7V H i > 19.9. 
The difference in abundance needed in order to exceed a 
given threshold of A^Mgii is thus not substantial between 
our DLA and sub-DLA samples. 

Fig. 2(a) plots [Zn/H] vs. Wtf$ for DLAs and sub- 
DLAs at z < 2 for which both quantities are available. 
(We consider only z < 2 here because Mg II A2796 is 
not available for the z > 2 systems.). Systems with large 
WJUg appear to show a wide range of [Zn/H], and sys- 
tems with high [Zn/H] also show a range of W^ral- ^ n 
other words, Mg II selection does not appear to bias the 
sub-DLA sample toward only solar or supersolar systems. 
Note that the fact that the sub-DLAs in Fig. 2(a) have 
[Zn/H] > — 1 merely states the observed difference be- 
tween DLAs and sub-DLAs that we have mentioned all 
along in this paper. It does not necessarily indicate a 
metallicity bias in the sub-DLA sample. 

To quantify the above considerations further, we now 
estimate what minimum [Mg/H] can be covered by the 
low-z sub-DLAs in our sample where the lowest W2796 
is 0.48 A. We use Mg II/ H I to estimate [Mg/H], given 
that (a) most photo-ionization calculations for sub-DLAs 
predict Mg/H to differ from Mg II /H I by at most 

0. 1 — 0.2 dex for the relevant range of ionization pa- 
rameters (e.g., Dessauges-Zavadsky et al. 2003; Meir- 
ing et al. 2006b); and (b) this also appears to be the 
case for diffuse interstellar clouds within the Milky Way 
(e.g. Jenkins, Savage, & Spitzer 1986; Cartledge et al. 
2006). If one assumes a linear curve of growth, then 
^2796 ^ O-^ A would correspond to log A Mg n > 13.1, 

1. e., [Mg/H] > -2.2 for log N m = 19.7. Of course, it is 
not correct to assume a linear curve of growth for such 
strong Mg II lines. To make more accurate estimates of 
the [Mg/H] threshold of our low-z sub-DLA sample, we 
generated a range of synthetic Mg II A 2796 Voigt pro- 
files with W2796=0.5 A. We considered 3 types of com- 
ponent velocity structures: (a) To begin, we considered 
the simplest and most conservative case of only 1 compo- 
nent (extremely unlikely given the complex Mg II A 2796 
line profiles observed in existing high-resolution stud- 
ies of DLAs/sub-DLAs). For a single component with 
an effective Doppler b parameter of 30 km s _1 , we get 
W2796 > 0.5 A if log N Mg n > 13.3, i.e. [Mg/H] > -2.0 
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for log TVhi = 19.7. For a single component with ef- 
fective b of 20 km s -1 , we get [Mg/H] > —1.7 for log 
TVhi = 19.7. (b) Next, we considered a more realistic, 
yet simple, case of multiple narrower components. For 3 
components with equal NMgii at velocities v = —50, 0, 50 
km s _1 with b = 10 km s _1 each, one would need total 
log N Mg n > 13.3, i.e., [Mg/H] > -2.0 for log N m = 19.7 
to have W^fgl > 0.5 A. (c) Finally, we produced profiles 
with 3 to 9 components with random relative A^gii, b, 
and v values. Each component was given a b value rang- 
ing from 5 to 30 km s" 1 chosen with a random number 
generator. The velocities of the components were taken 
to range from -100 to 100 km s _1 , again selected with a 
random number generator. Each component was given 
a starting column density distributed randomly in the 
range 9 < log Nm s ii < 11. Then, the column density of 
each component was increased by a small constant factor 
in successive iterations until the total EW of the profile 
was 0.5 A. Simulations were repeated for 3 to 9 compo- 
nents. Finally, the simulations were repeated for veloc- 
ity ranges of -200 to 200 km s _1 and -50 to 50 km s _1 . 
For each case, 1000 profiles were generated as described 
above, and histograms of the total Mg II column densities 
were examined. In these simulations, we obtained total 
median log A^gii in the range of 13.3 to 13.2 (with the 
total -/VMgii decreasing with increasing number of com- 
ponents). In other words, one can reach [Mg/H] > —2.0 
for W^l > 0.5 A for log N m = 19.7. Thus, for all ve- 
locity structures considered, the sub-DLA sample should 
include systems with [Mg/H] > —2, not only those with 
[Mg/H] > -1. 

Of course, systems with single narrow components 
would only be included in the W^gl > 0.5 A sample for 
higher metallcities, but in general Mg II A2796 lines in 
DLAs/sub-DLAs are not observed to show such narrow 
velocity structures. The Nm s ii needed to reach higher 
^2796 can b e larger, but such systems usually also have 
larger velocity spreads than systems with lower W^ral- 
For W^l = 1 A, our simulations with 5 to 9 compo- 
nents spread over -100 to 100 km s" 1 or -200 to 200 
km s _1 give [Mg/H] in the range of -1.5 to -1.7 for log 
N m = 19.7. For W£f e = 1.5 A, the simulations with 
5 to 9 components over the velocity range -100 to 100 
km s _1 give [Mg/H] of about -1.0 to -1.2 dcx for log 
N m = 19.7 or -1.3 to -1.5 dex for log N m = 20.0 (64 % 
of the sub-DLAs in our sample with W^g% > 1.5 A have 
log Am > 20.0). Even these values are almost an order of 
magnitude lower than the metallicities actually observed 
for most of these sub-DLAs. In any case, the sub-DLAs 
in our sample were not specifically chosen to have very 
large W^fgg- This is the resaon we calculate the [Mg/H] 
thresholds using the lowest W^g% m our sample (0.48 A). 
We therefore believe that while a bias against very low 
metallicity could be present in some situations (e.g., nar- 
row single components or very large W^fgl), this is not 
likely to explain the difference between DLAs and sub- 
DLAs in Fig. 1. Indeed, as Fig. 2(a) shows, the higher 
W2796 systems do not have systematically higher [Zn/H], 
and vice versa. Furthermore, as we show in section 2.3, 
our results would not be affected much even if there were 
a bias against sub-DLAs with [Zn/H] < —1. We note, 
however, that systematic observations of a much larger 
sample of DLAs and sub-DLAs are essential to study the 



issue of selection effects more definitively. 

2.2.2. Observational Sensitivity for Weak Zn II Lines 

Fig. 2(b) plots log iVznii vs - ^2796 f° r tne systems 
from Fig. 2(a). The DLAs show a narrower range of 
Aznii than the sub-DLAs. In fact, the DLAs and sub- 
DLAs with comparable W^fgg do not necessarily have 
similar Aznii- There are many ATznii detections or up- 
per limits for sub-DLAs that are actually below the DLA 
points at comparable WJy|g. Clearly, the sub-DLA mea- 
surements have adequate sensitivity to see low Zn II col- 
umn densities. The metallicities of some of the sub-DLA 
systems, derived after the fact by dividing the A ZnII by 
Ahi, are high (see Fig. 2a). The N m values play no di- 
rect role in the detection of Zn II lines. Of course much 
stronger Zn II lines could have been detected as well, and 
would in fact, have to be there if there were a significant 
number of solar or supersolar metallicity DLAs. How- 
ever, such strong Zn II lines for DLAs are not seen. On 
the other hand, there are a small fraction of sub-DLAs 
that have ATznii larger than those for many DLAs, and 
these are of course highly metal-rich systems. Thus, the 
larger fraction of metal-rich sub-DLAs does not appear 
to be a result of sampling similar iVznii values at lower 
Ahi values compared to the DLAs. 

2.3. The Role of Metal-rich and Metal-poor Systems 

It is interesting to compare the mean metallicities of 
the subsets of the DLA and sub-DLA samples after re- 
moving the systems with [Zn/H] < —1, since there can- 
not be a selection bias against the remaining systems. 
8 These subsamples consist of 50 DLAs and 27 sub- 
DLAs. Dividing these subsamples into 3 bins each, the 
iVni-weighted mean metallicity is —0.66 ± 0.15 in the 
lowest redshift DLA bin, while it is 0.08 ± 0.28 in the 
lowest redshift sub-DLA bin. While these numbers are 
certainly not representative of the overall observed DLA 
and sub-DLA populations (since by design they exclude 
the metal-poor systems), the DLA and sub-DLA means 
still seem to differ at 2 — 3 a level. Thus, the difference 
between DLA and sub-DLA metallicities at low z would 
not be affected much even if the sub-DLA samples were 
biased against systems with [Zn/H] < — 1. 

In principle, if there were a large number of metal-poor 
sub-DLAs with metallicity < 1/100 solar that could be 
missing from the Mg Il-selected samples, then the "true" 
sub-DLA mean metallicity could be lower. To assess this 
possibility, we repeated the survival analysis calculations 
after adding a hypothetical metal-poor sub-DLA popula- 
tion with a typical log Am = 20 in the lowest redshift bin 
(where the difference between DLA and sub-DLA mean 
metallicities is largest). We conclude that one would need 
to add 64 sub-DLAs with [Zn/H] = -2.0 (or 75 sub-DLAs 
with [Zn/H] = -1.5) in order to bring the Ani-weightcd 
mean metallicity of the sub-DLAs down to the DLA level 
(-0.82 dcx) in the lowest redshift bin. Such a large popu- 
lation of metal-poor sub-DLAs cannot be ruled out, but 
seems quite unlikely, given that there are currently only 
10 sub-DLAs with Zn measurements in this redshift bin 
(of which 6 are solar or supersolar). Turning the ar- 
gument around, if the fraction of metal-poor sub-DLAs 

8 Of course, given the discussion above, the sub-DLA samples 
should be sensitive well below [Zn/H]= —1 dex. 
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were the same as the fraction of DLAs with [Zn/H] < — 1, 
e.g., if there were 15 hypothetical sub-DLAs with a typ- 
ical [Zn/H] = -1.5 in addition to to the 10 observed sub- 
DLAs in the lowest redshift bin, then the TVni-weighted 
mean metallicity of the combined sample of 25 sub-DLAs 
would be -0.35 dex, not -0.82 dex as for DLAs. In other 
words, it would take a bimodal metallicity distribution 
very different from that for known DLAs to force the 
sub-DLA mean to equal the observed DLA mean. In 
fact, since the low-z DLAs are also discovered in the Mg 
Il-selected HST surveys, in principle there could be a 
population of low-z DLAs with [Mg/H] < —2.6 unde- 
tected in current DLA surveys that could drive the DLA 
mean metallicity to even lower values. Thus the DLA 
and sub-DLA means could still remain different. 

It appears that the primary reason for the difference 
between the observed DLAs and sub-DLAs in Fig. 1 is 
not a bias against metal-poor sub-DLAs, but the dearth of 
metal-rich DLAs, for which there should be no difficulty 
in detecting the metal lines. Fig. 1 shows that metal- 
rich systems can be found among the sub-DLAs, whereas 
they seem to be quite rare among the DLAs, despite 
much more extensive observations for the DLAs. Larger 
samples are needed, at any rate, to refine the differences 
between metal-rich sub DLAs and metal-poor sub-DLAs 
and DLAs, and those larger samples can be searched for 
more subtle or different selection effects than those we 
can study here. 

2.4. Other Considerations 

We note that owing to the relatively small size of the 
sub-DLA sample, the mean metallicities for this sample 
are sensitive to the particular details of the metal-rich 
sub-DLAs discovered so far. The maximum effect on 
the mean metallicity in the lowest- z bin comes from the 
sub-DLA with the highest N^aii in our sample, i.e., the 
system at z = 0.716 toward SDSSJ1323-0021 (Peroux 
et al. 2006a). We have adopted log N m = 20.21±°;^ 
for this system since this gives the best fit to the Ly- 
a absorption profile (Khare et al. 2004; Peroux et al. 
2006a). Raoet al. (2006) reported log N m = 20.54±0.15 
for this system, while Prochaska et al. (2006) reported 
log TVhi = 20.3. The value of Prochaska et al. is only 
stated to one decimal, without error bars, but is much 
closer to the value adopted here. In any case, if one were 
to take this system out of the sub-DLA sample, and put it 
instead in the DLA sample (adopting the Rao et al. 2006 
value for its Ahi), then the sub-DLA mean metallicity in 
the lowest redshift bin would become —0.26 ± 0.16. The 
corresponding value in the DLA sample would change 
only a little to —0.73 ± 0.15. (The DLA sample is larger 
and is less affected by an individual system.) Thus, the 
mean metallicity at low redshift would be still higher for 
the sub-DLA sample than for the DLA sample at the 3a 
level, but the difference would be smaller. 

Overall, we conclude that the higher mean sub-DLA 
metallicity does not appear to be an artifact of sample 
selection (on the basis of either Mg II or other metal 
lines), the use of Zn, or the use of Ani-weighting, and 
does not arise from insufficient Aznii sensitivity. How- 
ever, it is necessary to verify our results with (a) much 
larger sub-DLA and DLA Zn samples, and (b) more ac- 
curate Nm measurements for DLAs as well as sub-DLAs 
with higher resolution and higher S/N observations of 



the Ly-a profiles. 

3. CONTRIBUTION OF SUB-DLAS TO THE COSMIC 
METAL BUDGET 

Our knowledge of the cosmic metal budget is highly 
incomplete. Earlier estimates including contributions of 
high-z DLAs, Lyman-break galaxies, and Lyman-a forest 
suggested only ~ 10% of the amount of metals expected 
from the cosmic star formation history (e.g., Pcttini 
2004). Recent estimates also including sub- mm galax- 
ies and star-forming galaxies at z ~ 2 yield <~ 40 — 50% 
of the expected comoving metal density («4x 10 6 M Q 
Mpc -3 ; Ferrara, Scannapieco & Bergeron 2005; Bouche, 
Lehnert & Peroux 2005, 2006a; Pettini 2006). It has 
been suggested that dwarf galaxies contaminate the in- 
tergalactic medium with hot ionized gas (e.g., Tumlinson 
& Fang 2005), but the extent of this effect is not clear. 
Given the available observations, about half of the pre- 
dicted quantity of metals appears to be unaccounted for. 

In this context, it is interesting to quantify the contri- 
bution of the possibly metal-rich sub-DLAs to the metal 
budget. At low redshift it is difficult to estimate the sub- 
DLA contribution with the available data. The statistics 
for sub-DLAs and hence the comoving density of H I gas 
Ojji in sub-DLAs are not yet known at z < 1.7. However, 
one could assume that the relative H I contributions of 
DLAs and sub-DLAs at low z are similar to those at high 
z. At z ^2.5, the co-moving density of H I gas in DLAs 
and sub-DLAs is measured to be Oh i= 0.85 x 10~ 3 and 
0.18 x 10~ 3 , respectively (Peroux et al. 2005), i.e., a 
ratio O^p^/Ogy 4 = 0.21. Therefore, the co-moving 
mass density of metals in sub-DLAs, in units of O(Z ) 
(= ^baryons x Z = 5.5 x 10~ 4 ; Z Q =0.0126 by mass), 
is: n™ b - DLA = f x 10 xZ Q x O^- dla /O(Z )=/ x 
22.9 x J]™^ DLA where / (>1) accounts for the ionized 
fraction of the gas. For sub-DLAs near the high Ahi end, 
the ionized fraction is not important, but it may be im- 
portant for sub-DLAs with low JVhi. The exact value of 
/ will depend on detailed photoionization calculations as 
well as the Ahi distribution of the sub-DLAs, but is likely 
to be in the range 1-10. Thus, if the sub-DLA metallic- 
ity exceeds that of DLAs by a factor of about 7 at low z 
(Fig. 1), and the relative H I contributions of DLAs and 
sub-DLAs at low z are similar to those at high z, then 
the metal contribution of sub-DLAs will exceed that of 
DLAs by a factor of fx 1.5 [since Of LA = 2.6 x 10~ 3 , 
in units of O(Z )]. This factor may be even higher if the 
relative H I contribution of sub-DLAs at low z is higher 
than that at high z (e.g., the factor will be fx 4.5 if 
Ohi values of Prochaska et al. 2005 are used). These 
expectations will need to be confirmed with more data 
on sub-DLA Oh i and metallicity. 

At high redshifts, the mean global metallicity of the 
sub-DLAs is [(Zn/H)] « -0.7 (Fig. 1). Therefore, the 
co-moving mass density of metals in sub-DLAs, in units 
of O(Z ), is: n s z ub - D ' LA = f x 10-°- 7 x Z Q x 0.18 x 
lO- 3 /n(Z ) = f x 8.2 x 10~ 4 . Thus, at high redshifts, 
sub-DLAs contribute at least 32% of what DLAs con- 
tribute to the metal census. The sub-DLA contribution 
could be higher depending on the value of /. Prochaska 
et al. (2005) based on their analysis of the Sloan Digital 
Sky Survey (SDSS) Data Release 3 (DR3), find a value 
of Oh i for systems with 17.3 < logNni < 20.3 to be 
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0.57 times fin i in DLAs, a large fraction coming from 
the sub-DLAs. Thus, the sub-DLA contribution to the 
metal budget may be higher than that estimated above. 
We note, however, that the conclusions of Prochaska et 
al. (2005) will need to be confirmed at high resolution, 
which is essential to obtain reliable Nh i for sub-DLAs 
(Peroux et al. 2005). Taken together, the suggestions of 
high mctalicities at z <~ 2 for a few sub-DLAs and our 
finding of an even larger effect at low redshifts are consis- 
tent with sub-DLAs being more significant contributors 
to the total metal abundance of gas in the Universe than 
was previously realized. 

4. SUMMARY AND DISCUSSION 

While the sub-DLA samples are still small and not all 
the observed sub-DLAs are metal-rich, the fraction of ob- 
served metal-rich systems is larger for sub-DLAs than for 
DLAs. Our calculations suggest that the mean metallic- 
ity of the observed sub-DLAs is higher than that of the 
observed DLAs, reaching a near-solar level as expected in 
cosmic chemical evolution models. At z < 1, we find that 
the contribution of the observed sub-DLAs to the total 
metal budget may be several times higher than that of 
the observed DLAs. At higher z, the observed sub-DLAs 
may contribute > 32% of the DLA contribution to the to- 
tal metal budget. Furthermore, the results of Prochaska 
et al. (2005) on the relative number of sub-DLAs and 
DLAs may triple the contribution of sub-DLAs compared 
to the results based on Peroux et al. (2005). In general, 
our results support and extend the suggestions by Peroux 
et al. (2006a), York et al. (2006), and Prochaska et al. 
(2006) that sub-DLAs may contribute significantly to the 
overall global metallicity. 

One may wonder whether ionization effects are likely to 
be responsible for the difference we see between the global 
mctallicities for the DLA and sub-DLA samples. We be- 
lieve that this is not a significant effect for our samples. 
First, the Nm of the sub-DLAs in our sample are rela- 
tively high (86% systems with log N m > 19.7 and 73% 
with log Ahi > f9.9 in our primary sample.) Second, 
photoionization calculations presented in high-resolution 
studies of sub-DLAs (Dessauges-Zavadsky et al. 2003; 
Meiring et al. 2006b) show that while sub-DLAs have 
some ionized gas, the ionization corrections in abundance 
estimates are not significant, and would increase [Zn/H] 
if applied. Vladilo et al. (2001), Dessauges-Zavadsky 
et al. (2003), York et al. (2006), and Meiring et al. 
(2006b) all suggest that there is no trend in Aaiiii/Aaiii 
vs. Ahi in sub-DLAs and that the ratio is smaller than 1. 
Therefore, using Aznii/Ani in sub-DLAs for a measure 
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of Nzn/N-R does not introduce a systematic bias. 

It is important to establish whether the observed DLA 
and sub-DLA samples reflect their true properties or dif- 
ferential dust selection effects. If dust obscuration is sig- 
nificant and depends primarily on the metal column den- 
sity (e.g., Vladilo et al. 2006), then metal-rich DLAs may 
be more affected by dust selection than metal-rich sub- 
DLAs (Vladilo & Peroux 2005). Such an effect could 
make DLAs appear less metal-rich than sub-DLAs; but, 
this alone would not explain the large difference between 
DLAs and sub-DLAs seen in Fig. I, because so far there 
is no evidence for the necessary large number of high 
extinction DLAs (K06). K06 use this result and the 
observed mass-metalicity relationship for galaxies (e.g., 
Tremonti et al. 2004; Savaglio et al. 2005; Erb et al. 
2006) to conclude that sub-DLAs may arise in massive 
galaxies and DLAs in less massive galaxies, consistent 
with the relatively stronger evolution of metals in sub- 
DLAs compared to DLAs that we find here. That con- 
clusion is also consistent with trends in gas dispersion 
with metalicity of quasar absorbers (Lcdoux et al 2006, 
Meiring et al. 2006b). [But, see an opposing opinion 
on this last point by Bouche et al. (2006b), using a less 
direct argument.] 

The relatively smaller sample size for sub-DLAs com- 
pared to DLAs and the conflicting results on gas velocity 
dispersion just referred to emphasize the need for much 
larger samples of sub-DLAs and better velocity disper- 
sion measurements for the absorbers at a range of red- 
shifts, to confirm the results reached here. Large tele- 
scopes are essential to obtain adequate samples of excel- 
lent abundance and velocity dispersion measurements. 
The sub-DLA sample needs to be made at least as large 
as the DLA sample. It will also be important to obtain 
Ahi and [Zn/H] for absorbers toward UV-faint quasars, 
to confirm that there are not a large number of highly 
reddened quasars with foreground DLAs. More extensive 
observations with the HST Cosmic Origins Spectrograph 
to measure H I at z < 1.5 and Zn II at z < 0.55, and 
with ground-based spectrographs to measure Zn II at 
0.55 < z < 1.5 will be very important to confirm the 
findings presented here . 
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TABLE 1 

Zn Measurements in Sub-DLAs and DLAs 



Quasar 


Rcf. a 


z abs 


[Zn/H] 




Wznll/lO 1 * 


[X/H] 6 W 2796 


W26OO 










(cm -2 ) 


(cm- 2 ) 


(A) 


(A) 


Sub-DLAs: 
















0058+019 


1 


0.6125 


0.13 


1.20 + 0.42 


6.46 + 1.61 


1.63 


1.27 


1028-0100 


2 


0.6321 


-0.04 


0.79 + 0.26 


2.88 + 1.31 


1.58 


1.14 


1028-0100 


2 


0.7088 


< 0.06 


1.02 + 0.33 


< 4.71 


1.21 


0.89 


1323-0021 


3 


0.7160 


0.62 


1.62 + 0.78 


26.80 + 2.77 


2.23 


1.45 


0134+0051 


4 


0.8420 


< -0.36 


0.85 + 0.23 


< 1.48 


1.17 


0.86 



Note. — Tabic 1 is presented in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance 
regarding its form and content. 



TABLE 2 

Binned Sub-DLA Mean Metallicities for Various Samples 



Sample, 
No. of 
Systems 


Elements 
Used 


N m - 

weig- 
hted? 


Rcdshift Ranges of 3 bins" 


Mean Mctallicity in the 3 b 


ns 


l b , 30 
2 C , 32 
3 d , 30 
l e , 30 


Zn 
Zn 

Zn,S,Si 
Zn 


Y 
Y 
Y 

N 


0.61-1.14,1.15-1.87,1.89-3.17 
0.61-1.15,1.22-1.87,1.89-3.17 
0.61-1.14,1.15-1.87,1.89-3.17 
0.61-1.14,1.15-1.87,1.89-3.17 


0.05 + 0.27,-0.60 + 0.25,- 
0.02 + 0.25,-0.43 + 0.27,- 
0.05 + 0.27,-0.55 + 0.23,- 
0.16 + 0.19,-0.67 + 0.20,- 


-0.67 + 0.16 
-0.78 + 0.21 
-0.70 + 0.18 
-0.63 + 0.15 
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Fig. 1. — Logarithmic ^(H I)-wcighted mean Zn metallicity plotted vs. look-back time for DLAs and sub-DLAs. Dashed circles refer to 
the lowest time bin split into 2 bins with 10 DLAs each. The triangle denotes the formal lower limit to the average [Zn/H] for a composite 
spectrum from 698 absorbers with average log iV(HI) ~ 20 (sample 24) from York et al. (2006). The solid and long-dashed curves show, 
respectively, the mean metallicity in the chemical evolution models of Pei et al. (1999) and Somerville et al. (2001). 
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Fig. 2. — Left (a) [Zn/H] vs. rest equivalent width of Mg II A2796 for sub-DLAs and DLAs at z < 2 from our sample for which both 
quantities are available. Right (b): log N^nXl vs - ^2796 f° r * ne same se t of sub-DLAs and DLAs. (See text for further details.) 



